ON MULTIFRACTALITY AND TIME SUBORDINATION 
FOR CONTINUOUS FUNCTIONS 

STEPHANE SEURET 

Abstract. Let Z ; [0, 1] ^ K be a continuous function. We sliow tliat if Z is "lio- 
mogeneously multifractal" (in a sense we precisely define), tlien Z is tlie composition of a 
monofractal function g witli a time subordinator / (i.e. / is the integral of a positive Borel 
measure supported by [0, 1]). When the initial function Z is given, the monofractahty ex- 
ponent of the associated function g is uniquely determined. We study in details a classical 
example of multifractal functions Z, for which we exhibit the associated functions g and /. 
This provides new insights into the understanding of multifractal behaviors of functions. 



1. Introduction AND MOTIVATIONS 

Local regularity and multifractal analysis have become unavoidable issues in the past 
years. Indeed, physical phenomena exhibiting wild local regularity properties have been 
discovered in many contexts (turbulence flows, intensity of seismic waves, traffic analy- 
sis,..). From a mathematical viewpoint, the multifractal approach is also a fruitful source 
of interesting problems. Consequently, there is a strong need for a better theoretical under- 
standing of the so-called multifractal behaviors. In this article, we investigate the relations 
between multifractal properties and time subordination for continuous functions. 

The most common functions or processes used to model irregular phenomena are mono- 
fractal, in the sense that they exhibit the same local regularity at each point. Let us recall 
how the local regularity of a function is measured. 

Definition 1.1. Let Z e ^,"^^([0, 1]). For a > and to e [0, 1], Z is said to belong to 
Cj" if there are a polynomial P of degree less than [a] and a constant C such that, locally 
around tt^, 

(1.1) \Z{t)-P{t~t^)\<C\t~t^\^. 

The pointwise Holder exponent of Z at is hz{to) — sup{a > : / G }• 

The singularity spectrum of Z is then defined by dz{h) = dim{t : hz{t) — h] (dim 
stands for the Hausdorff dimension, and dim0 = — cx) by convention). 

Hence, a function Z : [0, 1] ^ R is said to be monofractal with exponent H > when 
hz{t) — H for every t £ [0, 1]. For monofractal functions Z, dz{H) — 1, while dz{h) = 
— oo for h H. Sample paths of Brownian motions or fractional Brownian motions are 
known to be almost surely monofractal with exponents less than 1 . For reasons that appear 
below, we focus on monofractal functions associated witli an exponent H e (0, 1]. 

More complex models had to be used and/or developed, for at least three reasons: the 
occurrence of intermittence phenomena (mainly in fluid mechanics), the presence of os- 
cillating patterns (for instance in image processing), or the presence of discontinuities (in 
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finance or telecommunications). Such models may have multifractal properties, in the 
sense that the support of their singularity spectrum is not reduced to a single point. Among 
these processes, whose local regularity varies badly from one point to another, let us men- 
tion Mandelbrot multiplicative cascades and their extensions ||6l [14] [l2l [T| , (generalized) 
multifractional Brownian motions ifTTl [3] and Levy processes (|4] [TOl (for discontinuous 
phenomena). 

Starting from a monofractal process as above in dimension 1, a simple and efficient 
way to get a more elaborate process is to compose it with a time subordinator, i.e. an 
increasing function or process. Mandelbrot, for instance, showed the pertinency of time 
subordination in the study of financial data 1.131 . From a theoretical viewpoint, it is also 
challenging to understand how the multifractal properties of a function are modified after 
a time change lfT9l ID . 

A natural question is to understand the differences between the multifractal processes 
above and compositions of monofractal functions with multifractal subordinators. 

Definition 1.2. A function Z : [0, 1] ^ K. is said to be the composition of a monofractal 
function with a time subordinator (CMT) when Z can be written as 

(1.2) ^ = .9°/, 

where g : [0, 1] ^ M is monofractal with exponent < < 1 and / : [0, 1] [0, 1] is an 
increasing homeomorphism of [0,1]. 

In this article, we prove that if a continuous function Z : [0, 1] ^ K has a "homoge- 
neous multifractal" behavior (in a sense we define just below), then Z is CMT. Hence, Z 
is the composition of a monofractal function with a time subordinator, and shall simply 
be viewed as a complication of a monofractal model. This yields a deeper insight into the 
understanding of multifractal behaviors of continuous functions, and gives a more impor- 
tant role to the multifractal analysis of positive Borel measures (which are derivatives of 
time subordinators). We explain in Section|6]and|2]how this decomposition can be used to 
compute the singularity spectrum of the function Z. 

Let us begin with two cases where a function Z is obviously CMT: 

1. If Z is the integral of any positive Borel measure [i, then Z = /c?[o.i] o Z, where the 
identity /d[o,i] is monofractal and Z is increasing. Remark that in this case, Z may even 
have exponents greater than 1 . 

2. Any monofractal function Zh can be written Zh = Zh o /(i[o,i], where Zh is 
monofractal and /(i[o,i] is undoubtably an homeomorphism of [0, 1]. 

These two simple cases will be met again below. 

To bring general answers to our problem and thus to exhibit another class of CMT 
functions , we develop an approach based on the oscillations of a function Z : [0, 1] — > R. 
For every subinterval / C [0, 1], consider the oscillations of order 1 of Z on / defined by 

UliZ) = sup \Z{t) - Z{t')\ ^ supZ(t) - inf 
t,t'G/ tei 

In tlie sequel, we assume tiiat Z is continuous and for every non-trivial subinterval / 

of [0, 1], uji{Z) > 0. This entails that Z is nowhere locally constant, which is a natural 
assumption for the results we are looking for 

It is very classical that the oscillations of order 1 characterize precisely the pointwise 
Holder exponents strictly less than 1 (see Section|2]i. 

Let us introduce the quantity that will be the basis of our construction. 
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Foreveryj > 1, fc e {0, 2-' — 1}, we consider the dyadic intervals /j^fe = [k2 (fc + 
1)2^-'), so that lJj,^Q 2J-1 -O.fe = ![' the union being disjoint. For every j > 1 and 
k e {0, 2^ - 1}, for simplicity we set cuj^kiZ) = lui. ^ ^Tj^l^) since Z is C"). 

Definition 1.3. For every j > 1, let Hj (Z) be the unique real number such that 

(1.3) 'j:mz))'/"^^'^^i. 

k=0 

We then define the intrinsic monofractal exponent of Z H{Z) as 

(1.4) H{Z)=\mi\TLiHj{Z). 

This quantity H{Z) characterizes the asymptotic maximal values of the oscillations of 
Z on the whole interval [0, 1]. This exponent is the core of our theorem, because it gives 
an upper limit to the maximal time distortions we are allowed to apply. 

It is satisfactory that H{Z) has a functional interpretation. Indeed, if Z can be decom- 
posed as (11.21 ). then the exponent of the monofractal function g shall not depend on the 
oscillation approach nor on the dyadic basis. In Section|4]we explain that 

(1.5) H{Z) = inf > : Z G S^^^^f ((0, 1))} = inf {p > : Z e ((0, 1))} , 

where -^^^'^^^^((0, 1)) and ©^''((0, 1)) are respectively the Besov space and oscillation 
space on the open interval (0, 1) (see Jaffard in ifTTI for instance). 

For multifractal functions Z satisfying some multifractal formalism, the exponent 
can also be read on the singularity spectrum of Z. Indeed (see Section HI, H{Z) corre- 
sponds to the inverse of the largest possible slope of a straight line going through and 
tangent to the singularity spectrum dz of Z. 

These remarks are important to have an idea a priori of the monofractal exponent of g 
in the decomposition Z — g o f. They also give an intrinsic formula for H{Z). 

Let us come back to the two simple examples above: 

1. For the integral Z of any positive measure /i, X)fc=o^ ^j.k{Z) — X)fc=o^ l^{^i,k) = 1. 
hence iJj (Z) = H{Z) — 1, which corresponds to the monofractal exponent of the identity 
/o?[o,i] from the oscillations viewpoint. 

2. The first difficulties arise for the monofractal functions Z^- When Zh is monofractal 
of exponent H, then we don't have necessarily H{Zh) = H. We always have H{Zh) < 1 
(see Lemma l23] in Section|2|i, but it is always possible to construct wild counter-examples. 
Nevertheless, we treat in details the examples of the Weierstrass functions and the sample 
paths of (fractional) Brownian motions in Section|5] for which the exponent H{Zh) meets 
our requirements. 

Unfortunately, the knowledge of H{Z) is not sufficient to get relevant results. For 
instance, consider a function Z that has two different monofractal behaviors on [0, 1/2) 
and [1 /2, 1]. Such an Z can be obtained as the continuous juxtaposition of two Weierstrass 
function with distinct exponents Hi < H2'- We have H{Z) = Hi, and Z can not be 
written as the composition of a monofractal function with a time subordinator This is a 
consequence of Lemma 12.41 which asserts that two monofractal functions gi and g2 of 
disctinct exponents Hi and H2 never verify gi = 52 o / for any continuous increasing 
function / : [0, 1] [0, 1] (indeed, such an / would "dilate" time everywhere, which is 
impossible). 
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We need to introduce a homogeneity condition CI to get rid of these annoying and 
artificial cases. This condition heuristically imposes that the oscillations of any restriction 
of Z to a subinterval of [0, 1] have the same asymptotic properties as the oscillations of Z 
on [0, 1]. 

Definition 1.4. Condition CI: 

Let J > 0, and e {0, V - 1}. Let Zj^k be the function 



where •^j,k is the canonical affine contraction which maps [0, 1] to Ij^k- 

Condition CI is satisfied for Z when there is a real number H > such that for every 

J > and X e {0, 2^ - 1}, H{Z,ik) = H{= H{Z)). 

Hence Zj^k is a renormalized version of the restriction of Z to the interval Ij^k- Re- 
mark that H{Zjji) does not depend on the normalization factor \/ljjj^k{Z). Although 
self-similar functions are good candidates to satisfy CI, a function Z fulfilling this condi- 
tion does not need at all to possess such a property. In order to guarantee that Z is CMT, 
we strengthen the convergence toward H{Zj^k)- 

Definition 1.5. Condition C2: 

Assume that Condition CI is fulfilled. There are two positive sequences (ej)j>o and 
{vj)j>o and two real numbers < a < f3 with the following property: 

(1) (£j),/>o and {rjj)j>o are positive non-increasing sequences that converge to zero. 



(2) For every J > and K e {0, ...,2'^ — 1}, the sequence {Hj{Zjj())j>i con- 
verges to H = H{Zj k) (it is not only a liminf, it is a limit) with the following 
convergence rate: For every j > [Jrjj], 



Assuming that H{Zj k) is a limit is of course a constraint, but not limiting in practice, 
since this condition holds for most of the interesting functions or (almost surely) for most 
of the sample paths of processes. Similarly, the decreasing behavior ilJi is not very 
restrictive: such a behavior is somehow expected for a C function. 

The convergence speed ( IL6b is a more important constraint, but the convergence rate 
we impose on (£,/)j>o toward is extremely slow, and is realized in the most common 
cases, as shown below. 

Tlieorem 1.6. Let Z : [0,1] Rbe a continuous function. 
Assume that Z satisfies CI and C2. 

Then Z is CMT and the function g in ([1.2i is monofractal of exponent H{Z). 

Remark 1.7. Such a decomposition is of course not unique: If Z is CMT and w : [0, 1] — > 
[0, 1] is C°° and strictly increasing, then Z — {g o w) o [w^^ o /), where 5 o it; is stiU a 
monofractal function of exponent H{Z) < 1 and o / is an increasing function. 

Nevertheless, if two decompositions il.2\ exist respectively with functions gi, 52, /i 
and /2, then gi and 32 are necessarily monofractal with the same exponent H{Z). This is 
again a consequence of Lemma lZ4l 





(1.6) 
(1.7) 



and for every fc e {0, 2^ — 1}, 



\H -H,{Zj,K)\<ej, 

< Lu,^k{Zj,K) < 2-^". 
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An important consequence of Theorem [L6| is that the (possibly) multifractal behavior 
of Z is contained in the multifractal behavior of /. More precisely, since / is an increasing 
continuous function from [0, 1] to [0, 1], / is the integral of a positive measure, say /i, on 
[0, 1]. The local regularity of /i is classically quantified through a local dimension exponent 
defined for every f e [0, 1] by 

= liminf H^^iifiM^ . I log. Mj^l^-^))! 

where B{t, r) stands for the ball (here an interval) with center t and radius r, and \ A\ is the 
diameter of the set A{\B{t,r)\ — 2r). The singularity spectrum of /i is then 

(1.8) d^{a) = dim{t : af^{t) = a\. 

It is very easy to see that if a^{ti^) — a, then hf{tQ) — aH. Hence for every h > Q, 
df{h) = d^{h/ H),\.t. there is a direct relationship between the singularity spectrum of Z 
and the one of /i. As a conclusion, Theorem I 1 . 61 increases the role of the multifractal anal- 
ysis of measures, since for the functions satisfying CI and C2, their multifractal behavior 
is ruled exclusively by the behavior of ^. 

As an application of Theorem I 1.61 we will prove the following Theorem 1 1.91 which re- 
lates the so-called self-similar functions Z introduced in ||9] with the self-similar measures 
naturally associated with the similitudes defining Z. 

Let us recall the definition of self-similar functions. Let be a Lipschitz function on 
[0, 1] (we suppose that the Lipschitz constant equals 1, without loss of generality), and 
let 5*0, 5*1, Sd-i be d contractive similitudes satisfying: 

(1) for every i ^ j, S'i((0, 1)) n ^^((0, 1)) = (open set condition), 

d-1 

(2) y ^^([0, 1]) = [0, 1] (the intervals Si{[0, 1]) form a covering of [0, 1]). 

We denote by < ro , r i , . . . , r^- 1 < 1 the ratios of the non trivial similitudes 5o , . . . , S'd- 1 • 

d-1 

By construction ^ = 1. Let Aq, Ai, Xd-i be d non-zero real numbers, which satisfy 



(1.9) < Xmin = min 

fc— 0,...,a— 1 



< Xmax = max 

fe=0,...,d-l 



< 1. 



Deflnition 1.8. A function Z : [0, 1] [0, 1] is called self-similar when Z satisfies the 
following functional equation 

d-1 

(1.10) ViG[0,l], Z(t) = ^Afc.(Zo(5fe)-i)(i) + 0(i). 

fc=0 

Relation ( |1.9t ensures that Z exists and is unique |'9l. 
Let us consider the unique exponent j3 > \ such that 

d-1 

(1.11) ^(A,)^ = l. 

This (3 is indeed greater than 1, since X]fe=o '^fe = 1 ™d l-^i I > for all k by ( 11.91 ). With 
the probability vector {jpo,Pi, ■■■,Pd-i) — (|Ao|'^, |Ai|^, |Ad-iP) and the similitudes 
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{Sk)k=o,...,d-i can be associated the unique self-similar probability measure satisfying 

d-l 

(1.12) fi = Y,Pk-{f^oS-'). 

Theorem 1.9. Let Z be defined by U.lOh Then, either Z is a n-Lipschitzfiinctionfor some 
constant k > (expliciteley found in Section^, or Z is CMT and there is a monofractal 
fiinction g of exponent 1/(3 such that 

(1.13) for every te [0,1], Z{t) = g{fi[0,t]), 

where ji is the self-similar measure ( I-/.72I ) naturally associated with the parameters used 
to define Z. 

The multifractal analysis of Z follows from the multifractal analysis of fj,, which is a 
very classical problem (see ISl). 

The paper is organized as follows. In Section [3] Theorem II. 61 is proved, by explicitly 
constructing the monofractal function g and the time subordinator /. Section |4] contains 
the possible extensions of Theorem 11.61 the explanation of the heuristics ( 11.51 ). and the 
discussion for exponents greater than 1 . In Section |5] |6] and |7] we detail several classes 
of examples to which Theorem 11.61 applies. First we prove that the usual monofractal 
functions Z with exponents H verify CI and C2. We prove Theorem 1 1.91 in Section |6l 
Finally we explicitly compute and plot the time subordinator and the monofractal function 
for a classical family of multifractal functions {Za)ae[o,i] which include Bourbaki's and 
Perkin's functions. 

Let us finish by the direct by-products and the possible extensions of this work: 
The reader can check that the proof below can be adapted to more general contexts: 

• the dyadic basis can be replaced by any 6-adic basis. 

• if (e,/) converges to zero (without any given convergence rate), then (under slight 
modifications of (?/,/)) the same result holds true. We focused on a simpler case, 

but in practice, a convergence rate e,j = o (^ (iogj)ii+« ^ shall always be always 
obtained. 

• The fact the the quantities H{Zj^k) are limits is only used at the beginning of the 
proof. In fact, only the existence of the scale [Jr]j\ such that ( 11.61 ) and ( 11.71 ) hold 
true at scale [Jrjj] is determinant. In particular, the conditions may be relaxed: We 
could treat the case where the H{Zj,k) are only liminf (and not limits). Again, in 
practice they are often limits, this is why we adopted this viewpoint. 

2. Preliminary results 

2.1. Oscillations and pointwise regularity. For every t e [0, 1], let Ij{t) be the unique 
dyadic interval of generation j that contains t, and I'^ [t) = Ij [t) + 2^^, IJ [t) = Ij {t) — 
2-3. 

Let us recall the characterization of the pointwise Holder exponents smaller than 1 in 
terms of oscillations of order 1 (see for instance Jaffard in IfTTI ). 

Lemma 2.1. Let Z : [0, 1] — > M a function, for some 7 > 0. Assume that hz{t) < 1. 
Then 
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In Lemma 1221 we impose some uniform behavior of the oscillations of Z on a nested 
sequence of coverings of [0, 1]. This is used later to prove the monofractality property 
of the function g in the decomposition Z ^ g o f (Section [3]), and also to decompose 
self-similar functions (in Section|6]l. 

Lemma 2.2. Let Z : [0,1] ^ R be a continuous function, t € (0, 1) and H e (0, 1). 
Suppose that there exists an infinite sequence (Tn) of coverings of[0, 1] such that 

• each Tn is a finite sequence of disjoint non-trivial intervals of [0, 1], such that 
UTeT„ T = [0, 1], 

• lim„^+oo maxTeT„ 1^1 = 0, 

• each interval T in Tn is contained in a unique interval T' o/T„_i, 

• for every T e Tn and T C T' & Tn-i, we have |T'|i+^" < \T\ < \T'\, for some 
positive sequence (Zn) that converges to zero when n +oo. 

Then: 

(1) If there exists a positive sequence (k„)„>i such that for every T G Tn, ujt{Z) < 
\T\"-'^", then for every [0, 1], hz{t) > H. 

(2) If there exists a positive sequence (k„)„>i such that for every T £ Tn, ujt{Z) > 
|T|^+«", then for every t e [0, 1], hz{t) < H. 

Remark that in part (2) of this Lemma, the property needs to be satisfied only for a 
subsequence {nk)k>i of integers. 

Proof. Let t E (0, 1), and r > small enough. For every n > l,t belongs to one interval 
T G Tn, that we denote T„(f). Denote by the smallest integer n so that T„(t) C B{t, r). 
By construction, t e T„^_i(t) and |r„,,_i(i)| > r (since r„,,_i(t) B{t,r)). By the 
fourth property of the sequence (T„), we have 2r > |T„^(<)| > |r„,,_i(i)|i+^"'- > 

Let us start by part (2), which is very easy to get. We have ujB(t,r) {Z) > ujt„ (t) {Z) > 
|r„^(t)|«+''". > r(i+^"^)(-f^+'=".-). 

Applying Lemma IZTl and using that and k„,, go to zero when r goes to zero, we 
obtain hz{t) < H. 

We now focus on part (1), which is slightly more delicate. If B{t, r) C r„^_i(t), then 

wehavea;B(t,r)(^) < c^T„_ut)(^) < |T„,-i(<)|^-'=--^ < (2r)(«--".-^)/(i+^--). 

If B{t,r) (t T„^_i(t), then there is an integer p (which depends on r) such that 
B{t, r)\Tn^-i{t) is covered by one interval T E Tp and not covered by any interval of 
Tp+i. Using the same arguments as above, we get |r| < r)\T„^_i(t)|i/(^+^p+i) < 
^i/(i+Zp+i) (remark that \B{t,r)\Tn^-i{t)\ < r). 

Now we have 

'^B(t,r)(Z) < LLlT^^_^(t){Z) + UJB(t,r)\Tr,^-i(t){Z) 

Since k„^, Zn^, Kp and Zp converge to as r ^ 0, Lemma lZTI vields hz{t) > H. □ 

2.2. Two easy properties for the study of H{Z). Let us begin with an easy upper-bound 
for H{Z). 



Lemma 2.3. Let Z : [0, 1] — > R /^e a non-constant continuous function. Then H{Z) < 1. 
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Proof. We can assume without loss of generality that lo\q^i\{Z) = 1. Let j > 1. By 

construction, ^^^^^ > 1. In order to have ( fT3] l. we necessarily have < 1. 

Hence the result. □ 

Lemma 2.4. Let gi and §2 be two real monofractal functions on [0, 1] of disctinct expo- 
nents < Hi < H2 < 1. There is no continuous strictly increasing function f : [0, 1] — > 
[0, 1] such that 51 = (?2 ° /• 

Proof. Suppose that such a function / exists. Let e > 0. This function / is Lebesgue- 
almost everywhere differentiable. There is a set E of positive Lebesgue measure such that 
for every t e E, f'{t) > 0. Around such a t, we have f{t + h) - f{t) ^ f'{t)h + o{t). 
Consequently, since hg^ {f{t)) — H2, for every \h\ small enough we have 

1(52 ° mt +h)- (52 ° /)(<)! < \f{t +h)- fit)\"^-' < c\h\"^-'. 

This shows that hg^of{t) > H2. Using again that hg^{f{t)) = H2, there is a sequence 
{h'n)n>i converging to zero such that for every n > 1, \g2{f{t) + h'^) ~ 92{f{t))\ > 
l^nl^^^"- Choosing /i„ so that f(t + /i„) = f{t) + h'^, we see that 

1(52 o f ){t + K) - (52 o f ){i)\ > + K) - /(t)|«^+^ > C\K\"^+^. 

This holds for an infinite number of real numbers (ft,„) converging to zero. Hence hg^of{t) ~ 
H2, which contradicts hg-^ {t) = Hi. □ 

2.3. A functional interpretation of H{Z). Note first that the previous results hold in the 
case where a 6-adic basis, > 2, is used instead of the dyadic basis. In fact, there is a 
functional interpretation of the exponent H{Z), independent of any basis, provided by the 
Oscillation spaces of Jaffard [,11 ] and the Besov spaces. Let us recall their definition, that 
we adapt to our context of nowhere differentiable functions. 

Let Z he a C function on (0, 1), where C is the global homogeneous Holder space 
and 7 > 0. Since |i9J where the theoretical foundations of multifractal analysis of functions 
were given, a quantity classically considered when performing the multifractal analysis of 
Z is the scahng function ijz{p) = sup |s > : Z e B'^{p^^{{0, 

Later, in 1 1 1 1, Jaffard also proved the pertinency in multifractal analysis of his oscilla- 
tion spaces Op^^{{0, 1)), whose definitions are based on wavelet leaders (we do not need 
much more details here). He also considered the associated scaling function Cz{p) = 

sup|s>0:ZeC';/P((0,l))}. 

Finally, still in lITTl . Jaffard studied the spaces Vp^^((0, 1)), which are closely related 
to our exponent H{Z), defined as follows: Denote, for j > 1 and k € {0, 2^ — 1}, 
ilj,k{Z) = i^ik2-3-32-3,k2-3+32~j]{Z), and consider the associated scaling function (we 
assume hereafter that Z is nowhere differentiable, as in Theorem lL6l l 

= l + limi^f^"g-^^-"'^"^--^^»^ 
For p > fixed, it is obvious that there is a constant Cp > 1 such that 

2^-1 2^-1 2^-1 

1/Cp ^ {u^jAzW < E (%fe(^))' < Cp E i^^MzW, 

k=0 k=0 k=0 
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since cjj,fe(Z) < flj.kiZ) < Y,ie{-3.-2,...,2.3} As a consequence, z/^ (p) = 1+ 

liminfj^+oo Sfc=o^^("j.fc(^))'' _ Comparing the definition of H{Z) with this formula, 
we easily see that H{Z) is the unique positive real number such that vz{l/ H{Zy) = 1. 

The main point is that the three scaling functions r/z, Cz and i^z coincide as soon as 
p > 1 [Tl |, andr]z{l/H{Z)) = Cz{l/H{Z)) = 1. Using the property of the Besov 
domains, we have 

H{Z)-^ = inf {p > : Z G ^^^^.""((0, 1))} = inf {p > : Z G Ol'^{Q, 1))} . 

2.4. Precisions for functions satisfying a multifractal formalism. Consider the scaling 
function C,z{p) above. Then for any function Z having some global Holder regularity |,1 IJ. 
Z is said to obey the multifractal formalism for functions if its singularity spectrum is 
obtained as the Legendre transform of its scaling function, i.e. 

for every /i > 0, dz{h) = inf {ph - Cz{p) + 1) (e U {-oo}). 

In particular, since (^z{^l H{Z)) = 1, we always have dz{h) < h/H{Z) (by using p = 
1/H{Z) in the inequality above). 

Moreover, assume that he — C,'z{^/H{Z)) exists and that Z satisfies the multifractal 
formalism associated with C,z at the exponent he- This means that the inequality above 
holds true for h = he, i.e. dz{hc) = hc/H{Z). 

From the two last properties we get that 1/H{Z) is the slope of the tangent to the 
(concave hull of the) singularity spectrum of Z, as claimed in the introduction. 

3. Proof of the decomposition of Theorem I1.6I 

The functions g and F are constructed iteratively. First remark that since {rjj) con- 
verges to zero, one can also assume, by first replacing rjj by max(77j , 1/ log j) and then by 
imposing that {ijj) is non-increasing, that the sequence {rjj) satisfies: 

• for every j > l,rij > 1/logj, 

• ijVj) is now a non-decreasing sequence and jrjj +oo when j +oo, 

• (rjj) still satisfies ( fL6] l and ( fLTb . 

Assume that conditions CI and C2 are fulfilled. 

3.1. First step of the construction of g and /. The exponent iJ(Zo,o) = H{Z) — H is 
the limit of the sequence Hj{Z), so there exists a generation Jq > 1 such that for every 

J>Jo, \H~HjiZ)\<eo. 

We set Hq = Hj^^ (Z), and by construction we have 'Yl,k=o i^Ja-ki^)) ' = 1- 
We then define the first step of the construction of the function /: we set 

m^Y.^ioj„,k'{z)f'"'' + {u:j,Az)Y'"''i:2'''t-k) ifte/,„,fe. 

fe'=0 

This function /q is strictly increasing, continuous and affine on each dyadic interval. More- 
over, /o([0, 1]) = [0, 1]. Let us denote Uj^^^k the image of the interval /,/(,, fc by /o, for every 
k £ {0, 2-'° - 1}. The set of intervals {Uj^^k : fc e {0, 2"'« - 1} clearly forms a 
partition of [0, 1). One remarks that 

(3.1) Vfce{0,...,2^«-1}, \fo{Ij„^k)\ = \Uj,M^i^JoAZ))'^""- 
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The first step of the construction of g is then naturally achieved as follows: we set 

9o{y) = Z((/o)-i(y))forye [0,1], 
orequivalently .go(/o(i)) = Z{t) for t e [0,1]. 
This function go maps any interval Uj„^k to the interval Z{Ij„_k), and thus satisfies: 

cou,,A9o)^^JoAZ) = \Uj,M"°- 
As a last remark, there are two real numbers < a' < j3' such that for every k 
2-JaP'/Ho < \Uj„^k\ < 2-''°°''/^°. Without loss of generality, we can assume that 
a' — a and (3' = (3 (a and /? appear in condition C2) by changing a into min(Q;', a) 
and (3 = max(/3, /?'), so that 

(3.2) foreveryfc, 2--'''^/"° < \Uj„^k\ < 2--^°°'/"°. 

3.2. First iteration to get the second step of the construction of g and /. We perform 
the second step of the construction. Let us focus on one interval Ijg^K, on which we refine 
the behavior of /q. By condition C2 and especially ( 11.6b . we have 

(3-3) Hjon.ol'kiZj„j<)y^"' = 1, 

fc'=0,...,2'-'o''-'ol-l 

where Hi = H[jg^j^]{Zjg^K) satisfies \H - Hi \ < ej^. 

Let Ji = Jq + [JoVJo]' hence Ji — Jq — [jQijjg]. Remark that, by il.li . we have for 
every k' e {0, 2["''"'-^ol - 1} 

(3.4) 2-^'o^''o]P < \co[j,,^^^^^k'iZj„,K)\ < 2-[^«''-'ol". 

Now, remembering the definition of Zj^ k, we obtain that for every k' e {0, ...,2'^^~'^° — 
1}, 

LU K2-'i^ +k' {Z) 
(3-5) UJj^^jg^k'{Zjo,K) = T7^. • 

Consequently, ( 13.31 ) is equivalent to 

J2 {u^j,Az)Y'"^^{u^.i.j<{z)f'"\ 

fe=0,...,2''l-l:/j^,fcC/jo,K 

and thus 

fe=0,...,2Jl-l:/jj.fcC/.7o,K 

We now define the function fi as a refinement on /q on the dyadic interval Ijo.k- We 
set for every k e {K2-''~-'°, {K + l)2-^^-''° - 1} and for t G Ij^^k 

hit) = MK2-''°) 

+ J2 {-j.,Az)r"Ku^.i.,K{z)Y'"''^"''^ 

+ (c.j„fc(Z))i/«i iu;j,^K{Z)y/""-'^"' i2'H ~ k). 

This can be achieved simultaneously on every dyadic interval Ija,K, K e {0, 2''" — 
1}, by using the same generation Ji for the subdivision (indeed, condition C2 ensures that 
the convergence rate of Hj{Zjg^k) does not depend on k). The obtained function is again 
an increasing continuous function, affine on every dyadic interval of generation Ji . 
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Let us denote Uj-^^k the image of the interval Ij-^^k by /i, for every k E {0, ...,2^^^ ~ 1}. 
The set of intervals {Uj^ k '■ k G {0, 2'^^ — 1} again forms a partition of [0, 1). We get 

(3.6) Vfce {0,...,2^^ -1}, \Uj,M = {^Ji,k{Z)Y'"'{u^Jo,K{Z)f'""'^'"'- 

but the main point is that we did not change the size of the oscillations of /q on the dyadic 
intervals of generation Jo, i.e. fi{Ij„,K) = /o(^Jo.a')- 

The second step of the construction of g is realized by refining the behavior of go'- Set 

.gi(y) = Z((/i)-i(y))forye [0,1]. 

This function gi maps any interval Uj-^^k to the interval Z{Ij^_k), and thus satisfies: 

u:u,,A9i) = ^.hAZ) with \Uj,,k\ = {u;j,,k[Z)Y'"^{o^JiMZ)f'"''-""'. 
Finally, we want to compare the size of the interval Uj^,k with the size of its father 
interval (in the preceding generation) Uj„_k- For this, let us choose k E {0, 2'^i — 1} and 
K G {0, 2''°-l} are such that C Ij^.k (hence A: can be written /c K.2-'^^-''>+k' 
with k' e{Q,..., 2"'!-"'° - 1}). Then, by ( B31 ), 

\UjiM = {L^jMz))'^''"{L^ji-j„.k'iZj,^K)y^"^ 

= \Uj„k\i^Ji-J„M{Zj,,K)y^"'. 

Using ( 13.41 ) we get 

\UjiM > \Uj„,k\^-^''''''"^^/"' - \Uj,M2-^'^"'-'°^^/"K 

On the other side, we know by ( |12] | that | Uj„ , fe | < 2 " -^o a / Ho ^ hence 

ic/j„,..i>ic/,7,..i>ic/,7„..r+"''"^, 

where the left inequality simply comes from the fact that /j^ ,k C Ij„,k- 

3.3. General iterating construction of g and /. This procedure can be iterated. Assume 
that the sequences {Jp)p>i, {fp)p>i and {gp)p>i are constructed for every p < n, and that 
they satisfy: 

(1) for every 1 < p < n, Jp = Jp_i + [Jp-i'nj^_i] and \H - Hp\ < ej^^„ 

(2) for every 1 < p < n, fp is a continuous strictly increasing function, affine on each 
dyadic interval Ij^^k and if we set fp{Ij^,k) = C^jp.fc, then 

p-1 

(3.7) \fp[Ij^.k)\ = \Uj^.k\ = {u^j^.k{Z)f'"^ n {^.J^,K,Ak){Z)Y'"-~""-^\ 

m=0 

where Km{k) is the unique integer such that Ij^^k C Ij^,K„^{k)^ for m < p, 

(3) For every I < p < n, the set of intervals {Uj^^k ■ k E {0, 2'^'' — 1} forms a 
partition of [0, 1). 

(4) For every 1 < p < n, if C/j„,fc C Uj^_^^K^^i(k), then 

|C/j,_„K,_,(fe)|'+"-'-^^ < \Uj^^k\ < 

(5) for every 1 < p < n, for y G [0, 1], 5p(y) ^ Z{{fp)-^[y)) 

(6) foreveryl < p < n, for every A: G {0, 2^-1}, we have /„i(/c2-P) ^ fp{k2-P) 
for every p < m < n. 
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The last item ensures that once the value of fp at k2 p has been chosen, every fm, 
m > p, will take the same value at k2^P. 

To build fn+i and gn+i, the procedure is the same as above. We use Jn+i = Jn + 
[Jnilj„], and we focus on one interval Ij^^k- We have by ( 11.61 ) 

fc=0,...,2-^" + i"-^"-l 

where Hn+i ~ H^j^^^^j ^{Zj^^^) satisfies \H — i/„+i| < £j„. we have for every k' G 
{0,...,2["'"''''.l - 1} 

(3.8) 2-[^""'''.l'5 < |c.[j„,,j.fc'(^j„,K)| < 2-[-^"''-^"l". 

and 

.OQ^ ^7 ^ ^ J„+i, if 2-^^+1 -^"+fc'(^) 

The same manipulations as above yield 

n 

(3.10) 5: {u.j^,,MZ)Y'"^^^^\l{^JM^){Z)Y'"-"''--^ 

fc=0,...,2'^"+l-l:/j,^^j.fcC/j„,K ™=0 
n-1 

Then/„+i is a refinement on /„: For every fc e ^X2'^"+^-'^'^ , {K+l)2-'''+^- 
and fori £ Ij„+uk 

fn+l{t) = fniK2-'") 

k—1 n 

+ E (c.,„,,..,(Z))V«.>.. U(^-^M^')iZ))'/''--'^"--^ 

n 

+ {uj^^,MZ)Y'"--^ X{{ojj^^^^^^^){Z)Y'"--^/"-^-{2'-^H-k). 

Remark that for every (fc, A:') e {K2'^"+^~'^" , {K + l)2'^"+i-'^" - for every m G 

{0,....,n},K™(fc) =if™(fc')- 

This can be achieved simultaneously on every dyadic interval Ij^.k, K G {0, 2"^" — 
1}, by using the same generation J„+i for the subdivision. The obtained function is again 
an increasing continuous function which is affine on every dyadic interval of generation 

Jn+l- 

We then define g^+i by .g„+i(y) = Z((/„+i)"i(y)) for y e [0, 1]. Let Uj^^^^k be the 
image of the interval/j^^^j.fc by /„+i, for every fc e {0, ...,2''"+^ — 1}. This function g„+i 
maps any interval Uj^^-^^k to the interval Z{Ij^^-^,k), and thus satisfies; 

'^'^J„+i,<e(5n+l) = UJj^+,,k{Z) 

with 

n 

m=0 

At this point, all the items of the iteration are ensured, except the item (4). We prove it 
now. As above, let us choose fc G {0, 2'^"+i — 1} and K € {0, 2"'^ — 1} are such that 
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/j„+i,fc C /j„,if, and let k' e {0, 2"'"+!-''" - 1} be such that k = K.2-^"+^-''" + k' . 
We have by ([53 

(3.11) - |C/,/„,fc|(c^J„+i-J„,fc'(^j„,K))'/''"+^ 
Then, by ( I3l4] |. 

As above, since by ( 13.21 ) that |?7j„,fc| < 2~'^""/^", we have 

ic/j„,fei > ic/j„+,fci > if/,7„,fcr"'"'^"^. 

3.4. Convergence of {gn)n>o and {fn)n>o- The convergence of the sequence (/„) to a 
function / is almost immediate. Indeed, each /„ is an increasing function from [0, 1] to 
[0, 1], and by item (5) of the iteration procedure, for every J > 1, for every fc € {0,...,2^ — 
1}, fm{k2~^) is constant as soon as > j. 

Recall that for every m and k, |/m(^j„,fe)| — \U,j^^k\- By ( 13. lit , and using ( |3.4| l, we 
obtain that \Uj^^^^k\ < |f^j„,fc|2~("'"+i~'^")", and iteratively 

(3.12) for every m > 1, |[/,/„,fe| < C2"-'"", 

for some constant C. Hence the sequence {\Uj^^-^^k\)m>i converge exponentially fast to 
zero, with an upper bound independent of k. 
As a consequence, if m > n, then 

1 1 fn f m 1 1 oo 

< max \fm{Ijr^ k)\ < max \Uj^ k\ 
fee{o,...,2-'™-i} ■ fee{o,...,2-''"-i} 

< (72""''"". 

This Cauchy criterion immediately gives the uniform convergence of the function series 
(/„) to a continuous function /, whose value at each dyadic number is known as explained 
just above. The limit function / is also strictly increasing, since it is strictly increasing on 
the dyadic numbers. 

The convergence of the functions sequence (5„)„>o is then straightforward. Indeed, 
each /„ is an homeomorphism of [0, 1], and admits a continuous inverse f^^. We thus 
have, for every n > 1, g„ ^ Z o f^^. The series (/~^) also converges uniformly on 
[0, 1]. Since Z is uniformly continuous on [0, 1], converges uniformly to a continuous 
function g : [0, 1] [0, 1]. 

Remark that / also admits an inverse function /"^, and that g = Z o f^^. 

3.5. Properties of g and /. Obviously, / is a strictly increasing function from [0, 1] to 
[0, 1], which is what we were looking for. All we have to prove the monofractaUty property 
of g. This will follow from Lemma l272l 

It has been noticed before that if we set, for every n > 0, r„ = {Uj^^k ■ k e 
{0, 2"'" — 1}}, then every T„ forms a covering of [0, 1] constituted by pairwise dis- 
tinct intervals. We obviously have: 

• lim„^4.oo maxTeT„ l^j = (using the remarks of Section [34l above), 

• (r„) is a nested sequence of intervals, 

• by item (4) of the iteration procedure, if T e r„ C T' G T„_i, then we have 
|T'|^+^" < |r| < |T'|, , with Z„ = VJ^-i^^r-- This sequence (Z„) converges 
to zero, since [rjn) converges to zero and {Hn) converges to H. 
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In order to apply Lemma 12.21 and to get the monofractality property of g, it is thus 
enough to prove the last required properties, i.e. there is a positive sequence (k„) converg- 
ing to zero such that for every T e r„, |r|^+''" < LUrig) < \T\"-'^". 

For this, let n > 1 and T e T„. This interval T can be written Uj^ k for some k E 

{0,...,2-'"-l}.Wehave|C/,„,,H(..,„,,(Z))V^" Ull.V^Jr.,K^ik)iZ))'/''--'^"--^ 
by construction , and g{Uj^^k) — gn{Uj^.k) ~ ujj^.k{Z). We just have to verify that 
lUj^^kl^^""" < ujj^^^k{Z) < |J7j„,fc|^"'"", for some k„ > independent of k. 
We have 

1 1 1 

log|C/j„,fc| = — logcjj„,fc(Z) + ^(77 77 )logwj^^^^(fe)(Z) 



1 



Writing that ^^^^ 

< + ^ + o(ej„)), we obtain 



(3.13) 



log|C/j„,fc| 1 



loga;j„,fc(Z) if 



< 172 (ej„ + + o(ej,„)) < + o(e,/,„)) and 

/e obtain 

2 EI^=o(gJ^ + o(gJ^)) log^J^.jf^(fc)(-Z") 
^2 logwj„,fe(Z) 



< 



Let us denote dm.k = — log loj^ ,K^(k) (Z), and i/im — ej^ for every m and k. Comparing 
the last inequality with the desired result, all we have to show is that 

1 A\ y^rri— n '4'mdm.k „ , 

(3.14) "'^ " ' > when n +00, 

(*n,k 

indepently of k. 

This is obtained as follows: Start from Ji, that we suppose (without loss of generality) 
to be greater than 100. Recall that, by the remarks made at the beginning of Section[3] we 
assumed that for every n > 1, rjj^ > (log Jn)^^. Subsequently, every term J„ is greater 
than In, where {In) is the sequence defined recursively by In+i — ^n(l + l/logZ„) and 
li = 100. Let us study the growth rate of such a sequence. It is obvious that lim„^+oo in — 
+00. We set Vn = login- We have Vn+i = Vn + log(l + l/w„) > + (1 — £)/vn for 
every n, with e that can be taken less than 1 /4 since vi is large enough. In particular, since 
vi > V2, V2 > \/2 + (1 — e)/\/2 > \/3- Recursively, if we assume that w„ > \/n + 1, 
then Vn+i > \/n + (1 — e)/-Jn > \/n + 1. Hence the sequence (In) converges to +00 
faster than exp and thus faster than any polynomial n^. ( We could be more precise, 
and prove using the same arguments that the growth rate of (Z„) is exactly exp i/ri.) 

Let us now find an upper bound for Using the item (1) of condition C2, we have 

that Ej ~ o (log jyi+K^ ^ ■ Using the lower bound we found for J,„ with e chosen small 

enough, we get that V'm — o(n^^/^(^+'')) < 0(71^^^"/^). 

The crucial point is that J2m>i ^™ ^ +00. Now, since d^.k is a sequence increasing 
toward +00, rewrite the left term of ( 13.131 1 as X]m=o ^^-j^, where < < 1. By a 
classical Caesaro method, we get that (13.14b is true, independtly of k. 

This directly impUes, by dTlJl l. that independlty of fc, - 
some sequence k„ that converges to zero. 



We now apply Lemma lZ2l which impUes that g is monofractal with exponent H. 



< K„, for 
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4. Around Theorem I1.6I 

4.1. Possible extensions for exponents greater tlian 1. Let us finally say a few words 
about functions having regularity exponents greater than 1 . The presence of a polynomial 
in the definition dl.ll ) of the pointwise Holder exponent is a source of problems when an- 
alyzing the local regularity after time subordination. Indeed, suppose that a continuous 
function gi behaves like \t — tol" (0 < a < 1) around a point to, and that another contin- 
uous function g2 behaves like a{t — to) + \t — to\^^^ (a 7^ 0) around to = .91(^0)- Then 
hgi{to) — a, hg2{to) — 3/2, but hg^ogi{to) — ct, which is different than the expected 
regularity 3a/2. Applying the construction above and getting a decomposition of a func- 
tion Z as Z = g o f, because of such problems, we didn't find any way to guarantee the 
monofractality of g. 

This is related to the fact that, still for the just above toy example, oJB{to,r}{g2) ^ 
2ar when r is small enough, while one would expect i-^B{to,r){g2) ^ r^^^. The use of 
oscillations of order greater than 2 (so that t^^^^j^ r)(52) ~ r^^^) was not sufficient for us 
to prove Theorem ll.6l for exponents greater than 1. 

An unsatisfactory result is the following: If Z has all its pointwise Holder exponents 
less than M > 1, then W1/2M ° Z has all its exponents smaller than 1 {W1/2M is the 
Weierstrass function ( 15.1b monofractal with exponent 1/2M), and one shall try to apply 
Theorem ll.6l to this function. 

As a consequence, this problem is still open and of interest. 

5. The case of classical monofractal functions and processes 

It is satisfactory to check that classical monofractal functions verify the conditions of 
Theorem ll.6l and that the exponent H{Z) is actually equal to their monofractal exponent. 
The proofs below are also representative examples of the method used to get convergence 
rates for lo H{Z). 

5.1. Weierstrass-type functions. LetO <a< > a and 6 > 1 be three real numbers. 
Let w be a bounded function that belongs to the global Holder class C''((0, 1)). Consider 
the Weierstrass-type function 

00 

(5.1) Z(t) = 

fc=0 

By ISI, either the function Z is , or it is monofractal with exponent a. For w{t) — sin(t), 
we obtain the classical Weierstrass functions monofractal with exponent a. In fact, it is 
proved in |5J that, if Z ^ (which is our assumption from now on), then there is a 
constant C > 1 such that 

(5.2) C"i2-^" < cjj- fe(Z) < C2-J". 

As a direct consequence, C^^/" < Yl=a < C"^^"' and obviously H{Z) = 

a. 

Let us find the convergence rate of Hj (Z) toward H{Z). We are looking for a value of 

e > and for a scale Jq for which Yl^k=Q ^'^'hk{Z)Y^^'^^^^ > 1' for every j > Jo- Let 
j > 1. We have, by (ES, 

2^-1 
k=0 
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For e small, l/{a + e) = 1/a — e/a"^ + o(e), and thus our constraint is reached as soon 

as 1 < C-^/o'+e/a'+oie)2]e/a+oUei) ^ This leads to 

je/a + oije) + (log^ C){~l/a + ejc? + o(e)) > 0. 

There is a generation Jq such that the last inequality is realized by e = ^iSSa-^i. Subse- 
quently, one necessarily has l/Hj{Z) > l/{a + e) for every j > Jo, since 

k=0 k=0 

and the mapping h X]fe=o^(^j fe(^))^^'' increasing with h. Hence Hj{Z) < a + e. 

Using the same method, we obtain Hj [Z] > a — 2e for j > Jq. 

Finally, we have found Jq large enough so that for every j > Jq, \Hj{Z) — a\ < £q, 
where we have set = 2e = 41og2 C /{Jq). 

For every J > 1 and K G {0, 2^ — 1}, we easily get the same convergence rates 
of Hj{Zj^K) toward a from the self-affinity property of the Weierstrass functions. More 
precisely, fix J and K, and let j > J + I- Remark that by construction of Zj^k, we have 

ujj^j^k{Z.j,K) ~ ^^''^Jj l^lz)^^ • ^'"^ looking for a value of e for which 



k=0.,...,2'-l:Ij,kClj,K 



^j,k(Z) 



> 1, 



for every j large enough. By (15.2b (used two times), and remarking that there are 2^ 
dyadic intervals of generation j included in Ij.k, we get 

f ^J.kjZ) 2/(a+e)r,-(i-J)a/(a+£) 

^ \^jk{z) - ^ 

k=0.....2l-l:I,_kClj,K ^ ■ ^ 

The same computations as above yield that, if we impose 77.7 = l/logj J and ej = 
4(log2 C) (log2 J) / J, then for every j > J + Jrjj, Hj-j {Zj^k) <oi — ej. Similarly, we 
obtain Hj{Z)>a + 2ej for j > J + Jrjj. 

Finally, for every j > J + J7]j, \Hj-j{Zjji) - a| < e,/, and Sj = o(l/(log j)2+''). 

Consequently, the Weierstrass functions satisfy CI and C2 with H ~ a, and they are 
also monofractal from our viewpoint. 

5.2. Sample paths of Brownian motions and fractional Brownian motions. Classical 
estimations on the oscillations of sample paths of Brownian motions {Bt)t>o yield ([ 1 1 1) 

P(^,,fc(St) <j2-^/2^ < ^exp(-jV) 

{Bt)>-2-^/A < i^cxp(-jV8) 



Hence, by a classical Borel-Cantelli argument, with probability one, there is a generation 
Jc such that for every j > Jc, we have the bounds j2^^^^ < ujj^kiBt) < j2^^^^ for the 
oscillations. 
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The same computations as for the Weierstrass functions show that there is a generation 
Jo such that if j > Jq > Jc , then 

> 1 and 'Y^{co,M)y^^'/'-'-'o) < l, 



fc=0 fc=0 
Jo 



where ej^ > C^-^ (for some suitable constant C). As a consequence, \Hj [Bt) — 1/2| < 
£jo for e very j > Jq. 

The self-similarity property of Brownian motions yields that for every J > 1 and K G 
{0, ...,2^ - 1}, for every j > J/log J, \Hj{{Bt)jji) - 1/2| < ej, where sj = C'-^, 
for some constant C independent of J and K. We omit the details here, that can be easily 
checked by the reader 

Consequently, a sample path of Brownian motion satisfies with probability one CI and 
C2, Wi\hH{Bt) = 1/2. 

Similar estimations on the oscillations of fractional Brownian motions Bh of Hurst 
exponent h lead to the same almost sure result for the sample paths, which also satisfy 
almost surely CI and C2 with H{Bh) = h. 

6. Applications to self-similar functions: Theorem|1.9I 



We consider the class of self-similar functions defined in Definition 11.81 with the pa- 
rameters of Z and the contractions Sk satisfying ( 11.91 ). 

The multifractal analysis of such a function Z is performed in Here we are going 
to prove that, under the conditions ( 11.91 ) on the Afc and the Sk, Z is CMT, and that the 
multifractal behavior of Z can be directly deduced from this analysis. It is a case where 
our analysis provides a natural way to compute the singularity spectrum of Z. 

6.1. Preliminary results on the oscillations of Z. Let us introduce some notations: for 
every n > 1, for every (ei, £2, £n) G {0, 1, d — 1}", we denote Iei.e2,---.en the inter- 

o 

val S'gj o S'ej o ... o S^^ ([0, 1]). The integer n being given, the open intervals {Iei,E2,--.,en) 
are pairwise disjoint, and the union of the closed intervals Iei,e2, ■■■,€„ equals [0, 1]. Now 
fix an integer n > 1 and a sequence (ei,e2, ...,£«) S {0, 1, ...,d — 1}". The interval 
Iei.62,-- . has a length equal to r^-^r^^ ■ ■ ■ r^^ . Finally, by iterating n times formula dLlOK 
we get that for every t e Iei.62,---,e„, 

(6.1) Z{t) = A,, -A,, •••Ae„ •(Zo5-io5-i_^o...o5-i)(i) 

+ A,, • A,, . . . Ae„_, • o S-^^ o S-l^ o ... o S-'m 
+ ... 

+ A,, -A,, •(0o5-;o5-i)(t) 

+ A,, -(00 5-1) (t) 

+ 

Recall that Xmax is defined in il.9i . 

Proposition 6.1. Let n — ^^'^'^'^ ■ Then either Z is a n-Lipschitz function, or there is a 
constant C > 1 such that for every n > 1, for every {ei, S2, Sn) G {0, 1, ...,d — 1}", 

(6.2) • |A,, • A,, • • • A,J < {Z)<C- |A,, • A,, • • • A,J. 
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Proof. We first find an upper-bound for w/^ ^ (Z). We use the iterated formula ( 16.1b . 

Letnand (£1,62, ■•■,£«) G {0, 1, 1}". Remark that when t ranges in 

(5^^^ o S~^_^ o ... o S~^){t) ranges in [0, 1]. Hence the oscillation of the first term of ( 16.1b 
is upper-bounded by jAe^ A^^ • • • Ae„ | ■uj[a_i\{Z). 

Now, for every k G {l,...,n— 1}, when i ranges in /e^, 62, {^7^ ° "^^k-i ° •■• ° 
S~^^){t) ranges in /^^.^j ... g^. Using that is a Lipschitz function, we get that the oscil- 
lation of each term of the form X^^X^^ '"^sk ' ('/' ° '^7k ° '^ek-i ° ■■■ ° ^ei^)i't) is upper 
bounded by | A^^ A^j • •• A^^ Kr^^^^ • • -re^). Finally, we obtain using ( |1.9b 



W7, 



i.e, .AZ) < \X,,X,,-KJ + J2\^e,Xe,--Xe,\-{r,,^,-r,J 



A:=l 

n— 1 k 



fc=l J=l 

n-1 



1 + 



<Ci|A,,A,, •••AeJ, 



/c=l 



where Ci = 1 + ^^fe^J Xmax < +00. 

We now move to the lower bound. Assume that Z is not K-Lipschitz. There are two 
real numbers < to,tQ < 1 such that \Z{tQ) — Z{to)\ > (k + ri)\tQ — to\, for some 
T] > 0. Let n and (ei, £2, £n) G {0, 1, d — 1}". Let us call tn = Si o S2 o ■■■ o 
Sn{to) and 4 = o 5*2 o ... o Snivo)- We obviously have i„,4 e Iei,e2,...,e^, and 
thus — \Z[i'n) ~ Z{tn)\. Using again (16.1b . we get by the same lines of 

computations as above 

\z{0-z{t^)\ 

n-1 

> |Ae,Ae2---A,J • \Z{Q - Z{to)\ - J2 \Xe,X,,---X,J{r,,^,---r,,M ~ H 

1^0-^0 1 



/c=l 



> \X,,X,,■■■X,^■\Z{t'^)-Z{^)\ 

\z{Q-z{t^)\ 



1-1 k 



|A,|y|z(o-z(to)i 



> |AeiAe2---A 



1 - 



1 



K + ?7 



+00 

fc=i 



> C2- |A,,A, 



•AeJ, 



where C2 = \Z{i'^) - Z(<o)|(l - ^^^) > by assumption. 
Finally, ( 16.2b is proved with C — max(Ci, Cj"^). 



□ 



6.2. Comparaison of Z with a self-similar measure. In order to prove that the function 
Z (II. 10b satisfies our conditions CI and C2, we introduce a self-similar measure /i, whose 
multifractal behavior will be compared with the one of Z, and the notion of multifractal 
formalism. 

Let us consider the exponent /3 > 1 such that ( 11. lib holds and the associated self-similar 
measure /i defined by (11.12b /i = J2k=o Pk ' ° ^k^)- '^^^ ^^^^ where the similitudes 
do not overlap, it is easily checked that by construction, for every n and (ei, £2, £n) <= 

{0, 1}", Wehave/x(/ei,e2,...,e„) = Pe.Pe^- ■ ■ Pe„ = I Aej A^^ • • • A^^ p. 
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This class of measures has been extensively studied ll6l 171 [161 [Tsll . For instance, the 
multifractal analysis of /i is very well known. For this, let us introduce the so-called L'^- 
spectrum of fi defined by 

(6.3) ■.qeR^T^{q) = \im inf (j, q) , where ( j, q) = ^"g^Efc^o /(^fc)' _ 

J^ + OO —J 

We only recall the properties we need ir7l [T6l[T8l 



Proposition 6.2. (1) For every g g M, (q) is the unique real number satisfying the 

equation 'Y^'l^^^iPkY i^kY'^^''^ — 1- The mapping q T^{q) is analytic on its 
support. Moreover, the liminf used to define (g) is in fact a limit for every q such 
that T^{q) is finite. 

(2) There is an interval of exponents — [ctmim c^max] such that for every a G I^, 
df^{a) = (rp)*(a), where (Tp)*(a) := infggR(qa - rp(g)) is by definition the 
Legendre transform of t^,. 

(3) If a ^ If,, then {x : = a} = 0. 

(4) There is M > 1 such that for every j, k large enough, 2^^^^ < n{Ij,k) < 2^-'/*^. 

Part (2) above is known as the multifractal formalism for measures, when it holds. 

Let us come back to the function Z. The reader can check that such a function Z 
satisfies CI and C2, and is thus CMT. Here we propose a quick proof of Theorem 1 1.91 
especially adapted to this case. 

The aim is to prove that Z can be written Z = go Each dyadic interval /j ^ is included 

in one dyadic interval Isi,...,e„, and contains a dyadic interval /^^ e„.e„+i, such that 

^^<i Iei,...,e„.en+i Can be Written respectively /j' fc' and/jv, j,// withO < 
j' < C, for some constant C independent of j and k. Consequently, w/^^ ^ ^ (Z) < 
uji- !_ (Z) < uji^_^ (Z), and thus 

C-V(/e„...,s„.s„+J^/^ < u;,MZ) < Cfi{h,_,S^^. 

Using now the self-similarity properties of the measure and the open set condition, we 
see that maxfe(pfe) • n{Iei,...,e„..e„+i) > lJ'[Ij,k) and n{Ij.k) > miiifc(pfc) • ii{I^^,,,,^^J. 
Hence, combining this with the last double inequality, we obtain that for every j and k 

(6.4) c-^^i{IJ,kY'f' < uj,^k{z) < c^i{i,,kY'p 

for another constant C, i.e. Proposition 16.11 extends to all dyadic intervals. 

Let us check that Z satisfies conditions CI and C2. Remark that, because of (I6.4l i. 

2^ — 1 2^—1 2^—1 

k=0 k=0 k=0 

Let El > 0. Using part (4) of Proposition |6]2| to find upper- and lower-bounds for uuj^k{Z) 
uniformly in k, we obtain 

k=0 k=0 
2'-l 2^-1 

and i^^AZ)^^'' < E MZ)nCKhk)r > C^+-^2-^--^*^. 



fc=0 k=0 
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Hence, the same computations as in the case of Weierstrass functions lead to the following 
choice: for some Jq large enough, we set ei ~ ^^tj ' ^'^'^ '^^'^ ^^'^ every j > Jq, 
\H,{Z)-H{Z)\<e,. 

Let now J, K be two integers, e > 0, and focus on H{Zj^k)- The same computations 
as above and as in the Weierstrass case yield for j > J 

fe'=0 fc=0:/j,fcC/j,K 



(iCTt) ^ (cMltw^) > Then 



First notice that 
we remark that 

Combining these inequaUties we get 

k'=0 *:=0:/j,fcC/.7,A- ^ ' ' ^ 

Let us focus on Since /3 > 1, we have > f 47^) ^- This quantity 

is lower bounded by L^J"-^)^ for some constant L (uniformly in k and K), since the ratio 
of the /^-measures of a dyadic interval and its father (in the dyadic tree) is uniformly upper- 
and lower-bounded for our dyadic self-similar measure /i. Finally, we obtain 

El,, (7 \\l3-e ^ ^-/3-2e [^],k(^)) r(.i-J)£ 

{u,^.,^,,{Zj,K)f > C 

k'=0 k=0:Ij,kClj,K ' ' 

k-ar rr ^(^'^^^^ 

fc— U:ij_fc Cij,jf 

Hence if we fix j — 1/ logj J, then the sum above is greater than 1 as soon as j > 

J + [Jtjj] and e > 4/3 i°g g i°g2 J ^ jhus for j > J + [Jtjj], H,^j{Zj^k) - H{Zj,k) < 

1 ^ 1 I „,;tu ^ _ 8 log C log2 J 

—^<^+£.j With ej - j^iog^ ■ 

Similarly one shows that for j > J + [Jrjj], Hj^j{Zjj<[) — H{Zj^k) > ^ — £j, and 
C2 holds true for Z. 

Applying Theorem |1.6l vields that Z is CMT and can be written as Z = g o F, where g 
is monofractal of exponent 

6.3. Computation of the singularity spectrum of F. Applying directly the construction 
of Section[3] we find a function g monofractal with exponent and a strictly increasing 
function / such that Z = g o f. 

One can even enhance this result as follows. Following the proof of Section[3] we see 
that for every p > l,foreveryfce {0, 1, — 1}, 

< \f{Ij,,k)\<Klj,.k)'-''-, 
where (Kp)p>i is a positive sequence decreasing to zero and n is defined by ( |1.121 l. 
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Let us denote by F the integral of the self-similar measure /i, i.e. for t E [0, 1] F{t) = 
/i([0, x]). We claim that f = gi o F for some function g which belongs to C^^''([0, 1]), 
for every 77 > 0. Indeed, define for every t e [0, 1] gi{t) = f o F~^{t). This is possible 
since F is an homeomorphism of [0, 1]. 

By construction, for every p > 1, for every k e {0, 1, — 1}, gi{F(^Ij^^k)) — 
f o o F(lj^^k) = f{l.j^,k), thus by the inequahty above. 



where we used that [F{Ij^^k)\ — i^{Ij^.k)- Now the sets of intervals {Fi^j^^k)) '■ k £ 
{0, 1, 2 'p — 1}} obviously forms a covering of [0, 1] to which Lemma lZ2] can be applied 
with H ^ I. 

Finally, we find that Z = gogioF = g2oF, where g2 is clearly monofractal with 
exponent since g and gi are monofractal respectively with exponents 1//3 and 1 (the 
resulting function 32 is monofractal since the oscillations of g and gi are upper and most 
important lower bounded on every interval). 

7. An example of function satisfying C1-C2 in a triadic basis 

We recall the contruction of multifractal functions of [15], which somehow generalizes 
the Bourbaki's and Perkin's functions. 

Let us consider the function Za defined for < a < 1 as the limit of an iterated 
construction: Start from = t on [0, 1], and define Zl{t) recursively on [0, 1] by the 

following scheme: Suppose that Z^ is continuous and piecewiese affine on each triadic 
interval [k3^^ , [k + 1)3^^], fee {0, 3^ — 1}. Then Z^^'^^ is constructed as follows: On 
each triadic interval [fc3^^ , (fc + 1)3^^], Z^^^ is still a continuous function which is affine 
on each triadic subinterval [fc'3"(^+i\ (fc' + 1)^-^^+^^] included in [fc3"^ (fc + 1)3^^], 



This simple construction is better explained by the Figure]?] 

It is straightforward to see that the sequence (Z^)j>i converges uniformly to a contin- 
uous function Za as soon as < a < 1. Bourbaki's function is obtained when a = 2/3, 
while Perkin's function corresponds to a = 5/6. 

For a < 1/2, the function is simply the integral of a trinomial measure of parameters 
(a, 1 — 2a, a), hence its singularity spectrum is completely known. We are going to explain 
why the functions Za, when a > 1/2, satisfy our assumptions, and thus can be written 
as the composition of a monofractal function g (with an exponent H we are going to 
determine) with an increasing function. We will also deduce from this study the singularity 
spectrum of Za- 

For a > 1/2, the limit function Za is nowhere monotone. Let us compute the oscilla- 
tions of Za on each triadic interval. 



and 



Zi+\k3'^) 
Zi+Hik+l)3-n 



Zi{k3-^) 

Zi{k3-n + a[zi{{k + 1)3-^) - Zi{k3-^)) 
Zi{k3-') + (1 - a)(zi{{k + 1)3-^) - Zi{k3-^)) 
Zi{{k + 1)3-'^). 
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Figure 1 . Iterated construction of Za, from step j to step j + 1 



Remark first that tlie slope of on [0, 1/3] is 3a, it is — 3(2a — 1) on [1/3, 2/3] and 3a 
on [2/3, 1]. Iteratively, if j > 1 and k e {0, 3-' - 1}, we write k^-^ = J2i=i ^p3"^, 
with e {0, 1, 2}. Then the slope of on [^3"^ (A; + 1)3--'] is simply 

(3a)"'='^'°(-3(2a- l))"'='^'H3a)"''-^'^ = 3^(a)"'='^'''(-(2a - l))"'='^'i(a)"'='^-% 

where n^.j^i is the number of integers p e {1, ...,j} such that ^ i (for i = 0, 1, 2) in the 
triadic decomposition of k3^^. 

Let us consider the trinomial measure /ia of parameters ( , , ). Then it is 
obvious that the absolute value of the slope of Z^ on each triadic interval [k3~^ , (fc+l)3~-'] 
can be written as /ia([fc3^-', (fc + 1)3^^])3-' (4a — 1)-'. As a final remark, we also notice 
that the oscillations of Za on each triadic interval [kS"^ , {k + 1)3^-'] is the same as the 
oscillations of Z^ on each triadic interval [k3~^ , (k + 1)3^-'], which is equal to 3^-' times 
the slope, i.e. 

(7.1) A^,([fc3"-'',(fc + l)3-^])(4a-l)-''. 

Let g e K. Let us compute the sum of the oscillations of Za at generation j. We have 

3^-1 3^-1 

(7.2) (^[fc3-. ,(fc+l)3-.](^a))' = E (Ma([fc3--'', (fc + 1)3- ^])y3''^ ^-^sHa-l) ^ 
k=0 k=0 

Let us explain now how we easily compute the exponent Ha such that jlAj holds true. 
For a multinomial measure jia (in fact, for any positive Borel measure), it is very classical 
in multifractal analysis to introduce the functions and the scaling function t^^ (q) 

defined for q S M as (16.3b but in the triadic basis: 

. ^ V • f ( ^ U f ^ l0g3 Efclro'(/^a([fc3-^ (fc + 1)3-^]))^ 

T'Pal?) = limulf Tp„j((?), where T^^,j{q) = — ^ — . 

In our simple case, it is easy to see that for every j >1 and q e M 

= - log3(2(a)« + (2a - 1)') + glog3(4a - 1). 
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What matters to us is the value of q for which the sum in (17.21) equals 1 . Let us write 
this specific value (7 as 1/ H, for some H > 0. When this sum is 1, then we have 

3^-1 

3-.?V(l///)3,(log3(4a-l))/H ^ ^ {fiam-^ {k + 1)3-^] ) ) a^Oogs (4a-l))/H ^ ^ 

k=0 

Let Ha be the solution of the equation —t^{1/ Ha) + log3(4a — 1)/Ha — 0, which is 
equivalent to 

(7.3) 2(a)i/^° + (2a - if'"'^ = \. 

This solution is positive, unique, and strictly smaller than 1. Hence, in this case, the 
monofractal exponent Ha is defined through an implicit formula. 

In order to get the whole condition C2, it suffices to notice that any rescaled function 
{Za)j.K (as defined in ( IL6l l, but here with triadic intervals) is actually equal to Za (if Za is 
increasing on [fcS^-' , (fc+l)3^^]) orto Za(l — .) (if Za is decreasing on [fc3~-' , (fc+l)3~-']). 
Hence H{{Za)j,K) is a limit for every J, K, and is even constant equal to Ha- Thus C2 is 
satisfied. 

We can then apply Theorem 1 1.61 and Za is the composition of a monofractal function 
Qa of exponent Ha with an increasing function Za- For a = 2/3, we see that Ha = 
1/2 is the solution to ^3^, since 2(2/3)^ + (1/3)^ = 1. We have plotted in Figure |7] 
the Bourbaki's function /2/3, its corresponding time change F2/3 and the corresponding 
monofractal function 32/3 of exponent 1/2 such that /2/3 = 32/3 ° -^2/3- 

In this case, we can even go further and compute the singularity spectrum of Za- The 
trinomial measure satisfy the multifractal formalism for measures, i.e. the singularity spec- 
trum of fj, is given by the Legendre transform of r^^ : 

c^Mo(") = (T-MJ*(a) := inUqa- Tf,^{q)), 

forevery a e [- log3(a/(4a - 1)), - log3(2a - l)/(4a - 1)]. 

It is easy to see, using ( 17. U . that if fia has a local Holder exponent equal to a at a point 
to, then Za has at to a pointwise Holder exponent equal to a^/^" — log3(4a — 1). Hence 
the multifractal spectrum of Za is deduced from the one of [la by the formula 

AzM = 4„ ((h + log3(4a - 1))!/^) 

forevery/iG [(-log3(a))i/^»)-log3(4a-l),(-log3(2a-l))i/^")-log3(4a-l)]. A 
more explicit formula is obtained as follows: for every q e R, if a = r^^, then (a) = 

The singularity spectra of /2/3 and fi2/3 arc given in Figure|2] 

Finally, remark that the maximum of the spectrum is obtained for — iiT^^y{0)y^^"- — 
log3(4a — 1), and dz^ (eta) = 1- After computations, we find aa — ^(log3(a^(2a — 1)). 

Let us consider the value of uq such that aag — 1. Then a§(2ao — 1) = 1/27, i.e. 
54aQ — 27aQ — 1. When a > qq, the set of points t for which hz^ (t) > 1 is of Lebesgue 
measure 1, hence we recover the main result of fT5\: Za is differentiable on a set of 
Lebesgue measure 1. Here we obtain in addition the whole multifractal spectrum of Za- 

Acknowledgment. The author thanks Yanick Heurteaux for a discussion on the oscilla- 
tions of the self-similar functions. 




Figure 2. Top: Bourbaki's function /2/3 on the left, the muhifractal 
time change F2/3 in the middle, and on the right the monofractal function 
(72/3 of exponent 1/2 such that /2/3 = 52/3 ° ^2/3- Bottom: Singularity 
spectra of 112/3 on the left, of /2/3 on the right. 



References 

[1] Barral, J., Mandelbrot, B., Multifractal products of cylindrical pulses, Probab. Theory Relat. Fields 124(3), 
409^30 (2002) 

[2] J. Barral and S. Seuret, The singularity spectrum of Levy processes in multifractal time. Adv. Maths 14 (1), 
437-468, 2007. 

[3] A. Benassi, S. Jaffard, D. Roux, Elliptic Gaussian random processes. Rev. Mat. Iberoamericana, 13(1):19- 
90(1997). 

[4] J. Bertoin, Levy processes. Cambridge Univ. Press (1998). 

[5] T. Bousch, Y. Heurteaux, On oscillations of Weierstrass-type functions. Preprint. 

[6] G. Brown, G. Michon, J. Peyriere, On the multifractal analysis of measures, J. Stat. Phys., 66:3^, 775-790, 
1992. 

[7] R. Cawley, R. D. Mauldin, Multifractal decomposition of Moran fractals. Adv. Maths. 92:196-236, 1992. 
[8] U. Frisch, G.S. Parisi, Fully developed turbulence and intermittency, Proc. International Summer school 

Phys., Enrico Fermi, 84-88, North Holland, 1985. 
[9] S. Jaffard, Multifractal formalism for self-similar functions: Part I and II. SIAM J. Math. Anal. 28(4), 

944-997 (1997). 

[10] S. Jaffard, The multifractal nature of Levy processes, Probab. Theory Relat. Fields 114(2), (1999) 207-227. 
[1 1] S. Jaffard, Beyond Besov spaces: Part 1, J. Fourier Anal. Appl. , 10(3), 221-246, 2004. 
[12] Kahane, J.-P, Sur le chaos multiplicatif, Ann. Sci. Math. Quebec 9, 105-150 (1985) 

[13] B. Mandelbrot, A. Fischer , L. Calvet, A multifractal model of asset returns, Cowles Foundation Discussion 
Paper, #1164, (1997). 

[14] B. Mandelbrot, Intermittent turbulence in self-similar cascades: divergence of hight moments and dimension 
of the carrier, J. Fluid. Mech. 62, 331-358 (1974) 



MULTTFRACTALITY AND TIME SUBORDINATION 



25 



[15] H. Okamoto, A remark on continuous, nowhere differentiable functions. Proc. Japan. Acad. 81, Ser. A, 
2005. 

[16] L. Olsen, A multifractal formalism. Adv. Math. 116(1): 82-196, 1995. 

[17] R. Peltier, J. Levy Vehel. Multifractional Brownian Motion, Technical Report, 2645, INRIA (1995). 
[18] Y. Peres, B. Solomyak, Existence of simensions and entropy dimension for self-conformal measures. 
[19] R. Riedi, Multifractal processes. Long Range Dependence: Theory and AppUcations, eds. Doukhan, Op- 
penheim, Taqqu, (Birkhauser 2002), pp 625-715. 
E-mail address: seuret@univ-parisl2 . f r 

Laboratoired'Analyseet de Mathematiques Appliquees- UniversiteParis-Est- UFR Sci- 
ences ET Technologie - 61, avenue du General de Gaulle, 94010 Creteil Cedex, France 



